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We present a simple quantum mechanical model exploiting the optical potential 
approach for the description of collision damping in the reaction ir^He —> d'N near 
the threshold, which recently has been measured at TRIUMF. The influence of the 
open d'N —> NNN channel is taken into account. It leads to a suppression factor 
of about ten in the d' survival probability. Applications of the method to other 
reactions are outlined. 
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I. INTRODUCTION 

In this paper we consider inelastic final state interaction (FSI) between the neutron and 
the hypothetical d' resonance with the quantum numbers T = 0, = 0“ produced near 

threshold in the reaction ir^He —>■ d'N. This resonance has been suggested to explain the 
observed peculiar resonance-like behaviour of the pionic double charge exchange (DCX) 
on nuclei, -\- A ^ A' -\- , ranging from Li to '^^Nh at ~ bOMeV (see i.0i and 

references therein). Such a A^A^-decoupled dibaryon had been predicted by QCD-inspired 
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models In a sequence of recent papers this behaviour for the DCX reaction, on Ca 

isotopes was qualitatively reproduced by Nuseirat et al. [0 and, for the medium 

and heavy nuclei, in the generalized seniority model, by Gibbs and Wu [§. Still there is 
no conventional model, which could demonstrate to reproduce these low-energy phenomena 
consistently for all nuclei. An independent test of the idea would be investigation of the ppn~ 
invariant mass spectrum in the double pion production, pp ppTT~7r~^, since the quantum 
numbers of the resonance forbid strong coupling to np. Recent measurements at CELSIUS 
0 do not show an evidence of the resonance at the predicted level in contrast to data 


from ITEP taken at a higher proton energy [jl^. Thus the situation appears not yet fully 
settled. Another way of searching for the d' resonance is DCX on Helium isotopes, where 
d' should manifest itself as a threshold phenomenon. At energies around lOOMeU the total 
cross section of the reaction ‘^He{7r~^,7i~) shows an excess over conventional calculations 
which could be ascribed to the contribution of the d' resonance H. At the same time, the 
contribution of this resonance for both, He{TT~^, tt' -) 0 and ^He{7i , vr"'') [|T^, occurred 
almost an order of magnitude smaller than expected. The reason for the suppression of the 
production cross sections near threshold can be due to the collision damping d'N NNN, 
and in the present paper we study this hypothesis using a simple quantum mechanical model 
to describe the propagation of an unstable resonance in the nuclear medium. Although 
singularities near production threshold always attracted much attention and are described 
in many textbooks (see, e.g., IB 1^ and references therein), we are not aware of similar 
approaches to the hnal state interaction suggested before. 

The paper is organized as follows. In the second section we present a simple quantum 
mechanical model exploiting the optical potential approach to account for the effect of colli¬ 
sion damping. In the third section we estimate the spreading width of the d' in the nuclear 
medium, namely, the virtual pion exchange contribution to the spreading is calculated in 
detail. In the next, fourth, section we give our numerical estimates of collision damping 
and of the d' survival probability in the nuclear medium, and compare the latter to the 
experimental data. The hfth section contains our conclusions and outlook. 
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II. A SIMPLE QUANTUM MECHANICAL MODEL 


In this section we present a simple quantum mechanical model to describe the propagation 
of a resonance in the nuclear medium. As an example, we consider the behaviour of the 
hypothetical resonance d' which can contribute to the pionic DCX on ^He. 

Due to a very short interaction time of the pion with the Helium nucleus it is the sudden 
influence to be the appropriate mechanism for considering the tt^Hc d'N reaction [ITS 


Thus the d'N system is created at the initial moment, t = 0, with the wave function of 
Helium 'ipue and then separates with time as 


t) = '0(r - vf, 0) = - vt). 


( 1 ) 


V = 



with the velocity v dehned by the energy excess over the threshold e, 

where fi ~ ^rriN is the reduced mass of the d'N system. 

If the potential in the effective Schrodinger equation for the wave function ([^) contains 
an imaginary part, 

U{r) = Uo{T)+tUi{r), t/i(r) < 0, (2) 

then the probability to hnd the system in its initial state, d'N, decreases with time ap¬ 
proaching a hnite limit with the survival probability being simply 

\i>d'N{t 00 )|^ 


Aii.r-i) — 


Let us, for simplicity, consider the step-like form of the potential Hi(r), 


( 3 ) 


f/i(r) = 


Ui = const, |r| < Ri 


( 4 ) 


0, |r| > Ri, 

where the potential strength Ui can be related to the elastic s-wave d'N zero-angle amplitude 
using Born approximation: 

ImF(O) U,{r)d\ = —fiRlU,, (5) 

whereas, according to the optical theorem, 

ImF(O) = ~ (6) 

47r 47r 
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The problem of separation of the elastic and inelastic parts of the cross section in such kind 
of reactions is a very subtle question. One of the relevant worries concerns the possible 
broadening of the ground state, so that a slight change in the resonance position could have 
been erroneously arcsibed to the inelastic part. Bearing in mind that, in general case, such 
an effect could lead to an overestimate of the damping in Eq. (||), we still do not expect it to 
change our estimates dramatically. Indeed, from a posteriori estimate (see Eq. (|4^) below) 
we hnd the elastic part of the cross section to be an order of magnitude smaller than the 
corresponding inelastic part, and thus it appears beyond the accuracy of the present paper. 

From Eqs. d^), (I) one hnds 

jj 3 (Tin'^ 

^ Sn Rl ' 

Notice that we do not meet any constraints on the wave length of the resonance since, in 
Born approximation, the zero-angle amplitude (j^) does not depend on the energy of the d'. 

As a direct consequence of the Schrodinger equation one has the probability conservation 
law at any moment of time in the form: 

4 f p(r,f)dV+ f i{r,t)ds = 2 j f/i(r)p(r, (8) 

J J S •/V 

where 

p(r,t) = |'^/’d' 7 v(r,^)|^ j(r,t) = + C.C. (9) 

are the probability density and the probability flux through the corresponding surface S', 
respectively. 

The term on the r.h.s. of Eq. (|^) describes absorption of particles. In the absence of this 
term the density of particles would evolve in space and time, according to Eq. (® , as 

Po(r,t) = \iJd'N{r - vt,0)|^ = - vt)p = po(r - vt,0), (10) 


if quantum spreading of the initial wave packet is neglected. 

Eq. (H) becomes especially simple if the size of the wave packet, hereinafter called i? 2 , is 
small compared to the radius of the optical potential, Ri. Indeed, one can split the entire 
evolution time into two periods, 0 < t < T and t > T, where T = Ri/v is the moment of 
time when d' and N leave the region of interaction. Then the solution to Eq. (^ reads 


P(r,t) = 


po{r,t)e t<T 
Po(r,t)e“^^, t > T, 


( 11 ) 
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with 

r = 2 |f/i|. ( 12 ) 

Hence the survival probability is 

f2UiRi\ / 

n;_ = exp(^^J=exp(^-^ 

where we also used the relation (|^. 

The opposite limit, R 2 ^ Ri, leads to the expression similar to (|T3|) with Ri interchanged 
with R 2 : 

Wsurv = exp • (14) 

To see this let us consider the following anzatz for the probability density: 



p{r,t) = pt,{T,t)w(t), 


(15) 


SO that the survival probability dehned above is just Wgurv = w{oo). After integration over 
an inhnitely large volume V (bounded by an inhnitely remote surface S) in Eq. (|^) one 
hnds the following equation for the function w{t): 


2w{i) = -w(t) J 2|t/i(r)|pi,(r-vi, 0)o!V, m(0) = 1, 


(16) 


which can be easily solved with the result 


tc(f) = exp (^~ J j 2|f4i(i’)|po(r — vT, 0)(i^r^ 


(17) 


We further simplify our estimates and consider the step-like form for both functions, 
f/i(r) and po(r,t): 


tli(r) =-^0(^1 - |r|), Pq{yR) = yQ{R2-\Y-wt\), 1/1,2 = 2- ( 18 ) 

For the two limiting cases, -Ri 3> R 2 and Ri -C -R 2 , the integral / 2\Ui\pQ(Rr entering Eq. 
(1^) equals 

^ [ 0(Ri - |r|)0(i?2 - |r - vt'\)d\ = ^Q{Ru, - vt'), (19) 

rir 2 J Vbig 

where 

Rbig = max {i?i, R2} , Vbig = max {Hi, V2} , (20) 


and we have neglected the edge phenomena. 
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FIG. 1: The diagram for the virtual pion exchange contribution to the spreading width of the d' 
resonance in the nuclear medium. 

Now the integral over t' in Eq. (p!7D can be done trivially that gives for the survival 
probability: 

Wsurv = exp , (21) 


Eft 


big 


which coincides with Eqs. (0), (0). 

In general case, ~ the survival probability can be presented as 


Wsurv = exp 


^ iv.R 


y F(RiIR^)\, R=^Rl + Ki V = 


( 22 ) 


with F being a smooth function of i?i/i? 2 , E(0) = F{oo) = 1. 

III. VIRTUAL PION EXCHANGE CONTRIBUTION TO THE SPREADING 
WIDTH OF THE d' RESONANCE IN NUCLEAR MEDIUM 


In this section we estimate the spreading width of the d' resonance in the nuclear medium. 
Let us start from the virtual pion contribution described by the diagram in Fig. 1. 


As shown in Ref. |]^, at low energies there is only one Lorenz invariant structure 
describing the d'NNir vertex and, therefore, the amplitude of the decay d' —> NNtt can be 
written as 


M 


d'^NN-K — 


f 


2m M 


ltlC75(iT2f)u^7r, 


(23) 


where mi ,2 are bispinors and C = 7270 - The coupling constant / can be expressed through 
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the NNtt decay width T = Fpp^- + + F^p^o as 


/ 


USti^V^MT 


3rio{M — 2mN — m„ 



(24) 


where rriN, and M are the nucleon, pion and d' masses, respectively; rjo being the 
enhancement factor due to the NN FSI in this decay, 770 ~ 4 5. For F 0.5MeV (as 

deduced from the data on the DCX reactions to discrete levels i) / ^ 14. 

The invariant matrix element of the process d'iV NNN, when all initial and outgoing 
particles are on mass shell, can be written in the form 


= = ^5^h,C'75(r2r)h 


^ijk 


ijk 


4:771 N 


ijk 


{P - pkY - m 


-Uk'y^TU, 


(25) 


3 

i=l 

where g is the pseudoscalar NNti coupling, /4ti ^ 14.3, Sijk is the totally anti-symmetrical 
tensor, and indices i,j,k = 1,2,3 numerate the outgoing nucleons. As a result, the amplitude 
Af is totally antisymmetric under permutation of the hnal nucleons. See also Fig. 1 where 
the notations are explained. 

The matrix element contains three different contributions, hence, when squared, it 
produces three different diagonal terms and three cross terms. After integration over the 
three-particle phase space all diagonal terms equally contribute to the cross section. The 
same holds true for the cross terms. Hence, it is sufficient to consider only one diagonal 
term squared {e.g., {ijk} = {123}) and one of the cross terms {e.g., {ijk} x {i'j'k'} = 
{123} X {321}). 

The diagonal term {123} squared and averaged/summed over spins and isospins of the 
initial/hnal nucleons equals to: 


i.e. 

^in 


(|.^S23r) = 124^(piP 2 + m;^) 


Ugf 


(Pp3 - m%) 




(26) 


‘w [(P -p-if - 

Integration over the phase space is essentially simplihed if the initial particles are at rest, 
pfj. — (TTijy^ 0 , 0, 0) and = (M, 0,0,0), which means the limit n —»• 0. In this limit 
= ad'N^NNN oc 1/v, SO that the product remains constant. Thus, in this limit. 


(l-^S23r)=2444AZ(P3-m^) 


(fg? 


TTIn 


(M uinY + m% - 2(M -F mN)E 2 , 
{ml - 2m% + 2 mArP 3)2 


(27) 












i.e. (lA^SsP) 
mjv and 


depends only on the energy of one of the nucleons, which ranges between 


(M + mAf)^ — 3m% 
2{M + niN) 


(28) 


Let us consider now the FSI between the nucleons in the relative s-state. Since there 
are three identical nucleons in the exit channel, only two of them can be in the relative 
s-wave with the third one being in the relative p-wave. It is easy to see that the matrix 
element (nucleon number k being in the NNti vertex) is proportional to the three- 

momentum |pfc|, therefore this nucleon is in p-wave in the rest frame of the initial d' and 
N and, hence, in p-wave relative to the nucleons i and j in the d'NNir vertex (see Fig. 1). 
At the same time nucleons from the d'NNir vertex are in relative s-wave [0. It means 
that only FSI between the nucleons in the NNti vertex is important. It can be taken into 
account by multiplying by 'i/’q(0), where '0q(r) is the continuum NN wave function 

(q (pi — pj)/2) containing the s-wave NN scattering amplitude, if the d'NNir vertex is 
assumed to be point-like: 


M 


ijk 




1 + 


R 


-1 


-a, 1 - iqij 


(29) 


where R ~ 0.8/m |^, a* is the scattering length, and is the three-momentum of 
either nucleon, i or j, in their centre-of-mass frame. 


Qij 



[{M + niNY + 


2Ek{M + rriN)] - rn%. 


(30) 


The Coulomb effects, which are important for very small invariant masses in the pp- 
subsystem, can be neglected in the integrated cross section. 

With the FSI taken into account, the diagonal term has to be replaced by 




1 + 


R 


-1 




iQij 


(31) 


Finally, the summary contribution of all diagonal terms to the differential cross section 


IS 


dcidiag 


3 1 


(27r)^ 3! AttinMv 


{\Mi23\ [Q + P - Pi - P2 - Ps), 


( 32 ) 


where 3! accounts for the three identical particles in the hnal state. 
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As already mentioned, in the limit n —> 0 the integration over the phase space is simplihed 
since there is no angular dependence, and 

J + -P — Pi — P2 — Ps) = y {\Mi23\‘^)dE2dE3. (33) 

Now let us consider contributions of the cross terms {e.g., {123} x {321}): 

- (XSXS* = -3hf 




(34) 


X 


iP2P3)iPpi - m%) + ipiP2)iPp3 - m%) + m%{P Y^iPi) ~ {PiP‘i){PP 2 + ra%) - m 

[{P - psY - ml] [(P - pi)2 - ml] 

As u —> 0 this expression depends only on energies of the outgoing nucleons, 

- {A4S°2h4S* + Aigbt"!) = - 3 — 


m 


(35) 


X (Pi - ?7Pv)[P3(M + mjv) - m%] + (P 3 - m7v)[(P3 - mAr)[Pi(M + mjv) - m%] 


-M 


-(M^ - 4m%) - (M - Pi - P 3 )(M + miv) 


[(P - P 3)2 - ml] [(P - Pi)2 - ml] ’ 


and the FSI can be taken into account in the way described above, so that each has 

to be replaced by Afpfc, as in Eq. (p9|). 

The integration over the three-particle phase space is similar to the integration performed 
in Eq. (^3]), so that for the total inelastic cross section one hnds 


(Tin = — mb, (36) 

V 

where the contribution of the cross terms is positive and does not exceed 15% when u —>■ 0. 
The product ainV is a smooth function of the energy over the threshold, e, so that, e.g., for 
£ = 20MeV (which corresponds to the kinetic energy of the nucleon of about 30MeV in the 
d' rest frame) the contribution of the diagonal terms is only 10% larger than for u —> 0. 

To proceed further we consider the d' propagation in the nuclear medium, which we 
assume inhnite, for simplicity. Then, similarly to (H^), one has for the spreading width: 


Ts = 2|ImPA|- 


(37) 


Here Ua is the summary effective potential created by nucleons in the nuclear matter, 

Ni 

Ua = J2u,^PaV3{Uo + zU3), 

i=l 


( 38 ) 
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where pa = is the density of the nuclear matter, Ni = paVi being the number of 

nucleons in the interaction region of the volume Vi. Using Eqs. (0) and (|37D one easily 
hnds: 

= 2paVi\Ui\ = PA(JinV, (39) 

or, if the Pauli blocking factor, p ^ 2.2, is taken into account^, then 

Ts = pA{iyin/v)v. (40) 

With the help of the result (|36D one can estimate the spreading width to be about 7MeV. 
Notice, however, that the pion exchange (the diagram in Fig. 1) contributes only 30 -i- 40% 
of the observed spreading width in nuclei, he., there must be an additional mechanism for 
the reaction d'N NNN, the total spreading width being 


r, ^ 10 20MeU. 


(41) 


Such a mechanism could be cx-exchange as discussed in Ref. [|]. 

Now we can check how well the approximation of the inhnite nuclear medium works. 
To this end we have to ensure that the typical free-path length of the d' due to collision 
damping, ~ v/Tg, does not exceed the radius of the nucleus Ra {^tiR^pa = Na)- For 
DCX to discrete levels this condition is fulhlled since for P^ given by (|4l|) and v ~ 1/10, 
which corresponds to the energy £ 3MeV above the threshold, both lengths appear to be 
of the same order of magnitude, and L^' < Ra- 

Besides we can perform a posteriori check of the validity of the approximations made in 
Section II. Namely, we can justify Born approximation used in Eq. (1):^ 


pR^lUoA 


V mN^s _ 1 

dvr paRi 3’ 


(42) 


and the one used in Eq. (|^) when the elastic part of the total cross section was neglected. 
Indeed, for the scattering off the step-like potential (0) one has 


= —tip^RUUq + f//), ain 

y 



(43) 


^To estimate Pauli blocking in the reaction d'N NNN we simply limit the phase space for the final 

nucleons by the constraint pn > Pp, where Pp is the Fermi momentum. 

^Let us remind the reader that Born approximation is valid for |Uo,i| 1 /(m^i)' 




with the ratio 
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where we put \Uq 


rsj 


a el 2 2 d 3 ^ 0+^1 

ain 3 \Ui\ Qtt Pa 

\Ui\ and substituted the velocity u ~ 1/10. 


1 

10 ’ 


(44) 


IV. NUMERICAL ESTIMATES 


In this section we present the results of numerical calculations using Eqs. 
the estimate ( 0 )- Thus for the survival probability one, hnally, has 


"^surv 


exp 


/ vRTs\ 
V PaVv) 




0, O, and 

(46) 


where we put the function F{Ri/R 2 ) equal to unity everywhere for the sake of simplicity. 

In Fig. 2 we depict the survival probability (^5|) as function of the excess energy of the 
initial pion over the threshold for d' production using for the spreading width T^ = 20MeV 
and for the radius R the values 1.4/m and 1.6/m, respectively. For the excess energy of 
3MeU above the threshold we find that only about 10% of all created resonances survive, so 
that the suppression factor is of order ten. Finally, in Fig. 3, we give the “bare” total cross 
section of the d' production in the reaction n^He —>■ d'N as well as the same cross section 
multiplied by the factor of the survival probability. We hnd our theoretical predictions to 
comply reasonably well with the experimental data given in Ref. (dots with error bars 
in Fig. 3), where already the effects of collision damping had been discussed briefly. 


V. CONCLUSIONS AND OUTLOOK 

Using an optical potential approach we have shown that collision damping strongly de¬ 
creases the survival probability of d' in presence of a nuclear medium. Fven for ^He it reduces 
the d' production cross section near threshold by an order of magnitude, thereby leading to 
a reasonable agreement with the data. Another way to approach the problem of collision 
damping of the d' resonance near its production threshold would be a full coupled-channel 
treatment. This requires the consideration of (at least) three channels: d'N and NNN 

(one could even neglect the vacuum width of the d' considering it as a stable particle). How¬ 
ever, solving the coupled-channel problem requires a rather accurate knowledge of many 
reaction amplitudes (with = 1/2“) involved in this problem. In principle, this would 
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FIG. 2; Survival probability Wsurv versus the excess of the pion kinetic energy over the threshold 
(in MeV) for F^ = 20MeV and R = 1.6/m (the curve 1) and R = 1.4/m (the curve 2). 

allow us to find the correction to the d' production amplitude by summing up all possible 
contributions to the amplitude TT^He —>■ d'N including the route ir^He NNN ^ d'N. 

Let us note that such effects should equally renormalize the amplitude of the d' formation 
in heavier nuclei. However, it is just this amplitude (deduced from DCX transitions to 
discrete levels [jl|) that was taken as an input to predict the d' production off ^He 
The difference between these phenomena is, therefore, due to different details of the d' 
propagation following its initial production by the incoming pion. In this respect it appears 
noteworthy to compare with the situation of A and E production in the reaction pp —>■ KNY 
near threshold, where the observed [0 surprisingly small S production cross section is 
interpreted as being due to pp —>• K^pTP K^pK, he., due to a strong FSI between p and 
which transfers E° immediately into the energetically much more likely A — a situation 
very similar to that of d'N NNN discussed above for the DCX on Helium isotopes. 
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FIG. 3; The cross section for the d! production in the reaction ir^He —> d'N (in fibarn) versus the 
kinetic energy of the pion (in MeV). The curve 1 gives the “bare” result, without the damping 
effect taken into account |17]; the curves 2 and 3 represent the same cross section multiplied by 


the survival probability, Eq. (^), for Tg = 20MeV and R taking the values 1.6/m (the curve 2) 
and 1.4/m (the curve 3). Experimental data |jl^ are shown by the dots with error bars. 
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